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Abstract 

The  trajectory,  travel  time  and  relative  approach  velocity  of  a  pursuer  tracking 
a  prey  along  a  simple  curve  of  pursuit  are  deduced  using  basic  principles 
of  two-dimensional  kinematics.  While  such  curves  are  well  known  in  the 
mathematics  literature,  little  attention  has  been  paid  to  this  problem  by  the 
physics  educational  community,  despite  the  fact  that  it  has  abundant  physical 
applications.  It  also  makes  an  interesting  alternative  to  the  traditional  problems 
of  introductory  kinematics. 

(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

Imagine  that  you  are  the  commander  of  a  submarine  and  spot  an  enemy  warship  travelling 
perpendicular  to  your  line  of  sight  at  distance  H  with  constant  velocity  of  magnitude  V.  You 
fire  a  torpedo  having  constant  speed  v  at  the  ship.  The  torpedo  tracks  the  target,  so  that  it 
always  travels  directly  toward  the  ship.  What  is  the  trajectory  of  the  torpedo  through  the 
water?  If  the  torpedo  impacts  the  warship,  how  long  does  it  take  to  reach  the  target?  If  it 
does  not  reach  it,  what  is  the  closest  distance  that  it  gets  to  it?  What  is  the  relative  approach 
velocity  of  the  ship  along  the  torpedo’s  line  of  sight? 

This  is  an  example  of  what  is  known  as  a  chase  problem  [1]  and  the  resulting  trajectories 
are  called  curves  of  pursuit  [2,  3],  Problems  of  this  general  sort  are  of  interest  to  the  military 
community  [4]  and  to  video  game  designers.  Historically,  the  particular  problem  in  this  paper 
was  first  solved  by  the  French  mathematician  Pierre  Bouguer  in  1732,  and  may  have  been 
originally  proposed  by  Leonardo  da  Vinci  in  the  form  of  a  cat  and  mouse  chase  [5],  It  is  treated 
in  some  second-year  undergraduate-level  differential  equations  texts  [6],  However,  such  books 
typically  use  the  formula  for  arclength  to  obtain  a  differential  equation  in  the  first  and  second 
spatial  derivatives  of  the  coordinates  of  the  curve  and  then  make  a  change  of  variable  to  solve 
it.  The  resulting  nomenclature  and  analysis  do  not  tie  in  well  with  undergraduate  physics 
courses.  The  present  paper  instead  uses  time  derivatives  and  vector  components  to  recast  the 
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Figure  1.  Trajectories  of  the  torpedo  and  warship.  Point  A  is  the  location  of  the  torpedo  when 
fired;  point  B  indicates  the  initial  position  of  the  target;  point  C  has  coordinates  (jc,  ;y)  and  denotes 
the  current  location  of  the  torpedo;  and  point  D  with  coordinates  ( Vt ,  H )  is  the  current  location  of 
the  ship. 


solution  in  terms  of  the  more  familiar  concepts  of  velocity  and  acceleration.  It  then  becomes 
an  application  of  2D  kinematics. 

2.  First  integral  of  the  differential  equation 

Set  up  a  coordinate  system  with  the  torpedo  starting  at  the  origin  and  travelling  initially  in  the 
y-direction  and  the  warship  steaming  in  the  v-direction.  At  some  arbitrary  instant  in  time  t 
(after  firing  the  torpedo  at  t  —  0),  the  positions  of  the  ship  and  torpedo  are  as  indicated  in 
figure  1.  The  dashed  right  triangle  has  abase  of  length  Vt  —  x,  height  H  —  y  and  hypotenuse 
(representing  the  distance  to  the  target)  of  length  L.  The  ratio  of  the  rise  and  run  of  this  triangle 
is  equal  to  the  ratio  of  the  y  and  x  components  of  the  torpedo’s  velocity,  because  the  torpedo 
is  headed  directly  toward  the  ship,  so  that 

Vt  —  x  =  (H  —  y)  — .  (1) 

Uy 

Differentiating  both  sides  with  respect  to  time  and  rearranging  leads  to 


using  the  fact  that  v2  +  v2  =  v2  is  a  constant  and  introducing  the  y-component  of  the  torpedo’s 
centripetal  acceleration.  In  order  to  integrate  this  equation,  one  can  substitute  the  chain-rule 
result,  ay  =  uy  duy/dv,  familiar  from  the  derivation  of  the  work-kinetic-energy  theorem.  The 
variables  can  then  be  separated  and  both  sides  integrated  to  obtain 

v  r _&i_  =  _  n  vdvy  (3) 

v  Jo  H  —  y  Jv  UyJv2  —  Uy 


where  the  lower  limits  are  the  initial  values.  Computing  these  standard  integrals,  one  gets 
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The  quantity  in  the  last  set  of  parentheses  appears  so  frequently  in  the  remainder  of  this  paper 
that  it  is  useful  to  give  it  its  own  symbol,  z  =  1  —  y/H,  in  terms  of  which  equation  (4)  can  be 
rewritten  as 

2v 

Vy  =  z-v/v  +  zv/v  ■  & 

This  is  the  solution  for  the  v-componcnt  of  the  torpedo’s  velocity  along  its  curve  of  pursuit. 

3.  Solving  for  other  variables  of  interest 


Equation  (5)  can  be  used  to  readily  solve  for  other  physically  interesting  quantities.  They 
will  be  expressed  in  terms  of  the  dimensionless  parameter  z,  which  represents  the  normalized 
transverse  distance  from  the  ship  to  the  torpedo.  First,  v2  +  v2  —  v2  leads  to 


vx  1  -  z2v/v 

cos  6  =  —  =  - — . 

v  1  +  z2v 


(6) 


Equations  (5)  and  (6)  can  alternatively  be  compactly  expressed  in  terms  of  hyperbolic  functions. 

Next,  to  get  the  coordinates  of  the  torpedo,  note  that  vy  =  dy/dt  so  that  equation  (5)  can 
again  be  separated  and  integrated  to  get 

-V/v  /  ^  \  V/v~ 

+ 


(=ir 

2u  J o 


1- 


H 


1-*- 


H 


dy. 


( 1 ) 


Provided  that  u  ^  V,  the  integrands  are  of  the  form  zd 
t 


where  n  ^  1,  so  that  the  result  is 


y  +  ”  zi-v/«  +  ^ 


7\+V/v 


r  2v  2  u 

where  r  is  a  characteristic  time  defined  as 
Hv 

x  =  - 


+  1, 


(8) 


(9) 


v*  —  V 2 

whose  physical  significance  will  become  clear  momentarily.  Equation  (8)  gives  the  time  that 
the  torpedo  needs  to  reach  a  given  position.  One  can  substitute  that  along  with  equations  (5) 
and  (6)  into  (1)  to  find 

* _  0-5  Ti-v/v  +  0-5  „i+v/u _ V/v 

H~ 


(10) 


(11) 


V/v  -  r  V/v+l  ““  V2/v2  -  1 

which  is  the  trajectory  x(y)  of  the  torpedo  expressed  in  terms  of  the  normalized  quantities 
x/H,  y/H  —  1  —  z  and  V/v.  This  is  plotted  in  figure  2  for  several  values  of  the  speed  ratio. 
The  torpedo  strikes  the  target  when  y  —  H  ^  z  =  0,  which  implies 
HV/v 

■*imPact  _  i  _  v2/v2' 

This  is  positive  if  and  only  if  v  >  V.  As  expected,  the  torpedo  only  hits  if  it  is  travelling 
faster  than  the  ship !  In  that  case,  the  time  between  launch  and  impact  is  Ximpact/  V  =  x  defined 
by  equation  (9).  This  impact  time  r  can  be  reduced  if  the  torpedo  leads  [7]  or,  even  better, 
mirrors  [8]  the  ship’s  motion  rather  than  heading  straight  toward  it. 

The  distance  between  the  torpedo  and  target  at  any  time  is 

L  =  S[{VT 


s  V 

y)—, 


■x)2  +  (H  -  y)2  =  ( H 

uy 

using  equation  (1)  in  the  second  step.  Substituting  equation  (5)  then  gives 

7l+V/vy 


(12) 


L  =  f  (z1-^ 


+  z 


(13) 
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Figure  2.  Graphs  of  equation  (10)  or  (16)  for  0  y  ^  H .  The  number  on  each  curve  labels  the 
value  of  V/v.  Impact  occurs  if  V/v  <  1 ,  as  is  evident  for  the  top  two  traces. 


The  value  zmin  that  gives  the  distance  of  closest  approach  Lmm  for  V  >  u  is  found  by  setting 
dL/dz  =  0, 


v-vy/2V 

V  +  v) 


-  H  (V-v\v/2V 
yj\  —  v2/V2  V  V  +  V  J 


(14) 


As  one  might  have  guessed,  the  minimum  separation  is  equal  to  H  in  the  limit  as  V /u  -»  oo 
and  occurs  at  t  =  0. 

On  the  other  hand,  as  V  v,  equation  (14)  has  the  limiting  value  7,mm  —*■  H / 2  and 
this  occurs  as  t  —*■  oo.  Assuming  the  torpedo  does  not  actually  hit  the  target,  it  can  never 
get  closer  than  this,  regardless  of  its  speed !  The  trajectory  in  this  case  is  given  by  integrating 
equation  (7)  with  V/v  =  1  to  get 


t 


H 

4v 


(1  -z2-2lnz). 


Then,  following  the  same  reasoning  used  to  derive  equation  (10),  one  deduces  that 


x 

77 


In  sfz, 


(15) 


(16) 


plotted  as  the  middle  curve  in  figure  2. 


4.  Relative  approach  velocity  of  the  target 


From  figure  1,  one  can  see  that  the  ship  is  receding  from  the  torpedo  at  speed  V  cos  6  along 
the  line  of  sight,  while  the  torpedo  is  approaching  along  the  same  line  at  speed  v.  Hence,  the 
normalized  relative  approach  velocity  of  the  warship  is 


frel  .  V 

-  =  1 - COS  ( 


(17) 


where  cos  0  is  given  by  equation  (6).  This  is  plotted  in  figure  3  for  the  case  of  V  jv  =  2. 
A  positive  value  of  vrc\/v  indicates  that  the  ship  and  torpedo  are  approaching  one  another; 
a  negative  value  means  they  are  receding.  Therefore,  the  intercept  along  the  horizontal  axis 
represents  the  position  of  closest  approach  between  the  torpedo  and  target.  This  can  be  verified 
mathematically  by  substituting  /mln  from  equation  (14)  into  (6)  to  deduce  the  heading  of  the 
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Figure  3.  Plot  of  the  normalized  relative  approach  velocity  from  equations  (17)  and  (6)  as  a 
function  of  the  normalized  transverse  position  of  the  torpedo,  for  the  case  of  V  =  2v.  At  y  =  0, 
the  torpedo’s  heading  is  6  =  90°  and  thus  vrt\/v  =  1.  On  the  other  hand,  as  y  — >  H,  the  torpedo 
falls  in  line  behind  the  ship  (since  V  >  v)  so  that  6  =  90°  and  hence  vre\/v  =  —  1.  The  horizontal 
intercept  is  located  at  y/H  =  1  —  3-1/4  =  0.24. 


torpedo  at  minimum  separation, 


=  cos  — 
V 


(18) 


or  60°  in  the  case  of  figure  3.  This  simple  relation  can  alternatively  be  deduced  by  setting  the 
time  derivative  of  the  first  equality  in  equation  (12)  to  zero,  and  then  substituting  equation  (1) 
into  it. 

In  practice,  the  curve  in  figure  3  could  be  experimentally  mapped  out  if  the  torpedo 
measures  the  Doppler  shift  in  the  tracking  return  signal.  The  classical  Doppler  shifted  reflection 
frequency  off  the  target  is 

Xu  —  V  cos  0)  (u  +  v) 


f'  =  f- 


09) 


{u  —  v )  ( u  +  V  cos  6)  ’ 

for  a  source  of  frequency /aboard  the  torpedo  with  a  signal  speed  u  through  the  water.  (For 
example,  in  the  case  of  Doppler  sonar,/is  typically  of  the  order  of  hundreds  of  kilohertz  and 
the  sound  speed  u  is  about  1500  m  s-1  in  seawater.)  The  four  terms  in  parentheses  arise  from 
the  Doppler  shifts  in  the  signals:  emitted  by  the  moving  source — first  term  in  the  denominator; 
received  at  the  ship — first  term  in  the  numerator;  echoed  by  the  ship — second  term  in  the 
denominator;  and  received  back  at  the  torpedo — second  term  in  the  numerator.  For  realistic 
torpedo  and  warship  speeds,  u  and  V  are  much  smaller  than  u,  so  that  equation  (19)  can  be 
approximated  as 


/'  =  /  1 


2 —  cos  ( 
u 


„  u 
+  2  — 
u 


Af  u 
— -  =  2 
/  v 


1 


V 

—  cos  0 
u 


(20) 


where  A /  =/'  —  /.  The  left-hand  side  of  the  second  equality  is  the  normalized  frequency 
shift  divided  by  the  Mach  number  of  the  torpedo.  Meanwhile,  the  right-hand  side  is  twice 
the  normalized  relative  approach  velocity  of  the  target,  according  to  equation  (17).  Hence, 
figure  3  directly  gives  the  Doppler  shift  to  within  a  constant  scaling  factor.  In  particular,  zero 
shift  represents  the  best  opportunity  for  the  torpedo  to  explode  and  attempt  to  inflict  damage 
on  the  ship  (assuming  it  has  not  already  impacted  it  by  then),  which  is  the  basis  of  proximity 
fuses.  Similar  spectral  effects  are  exploited  by  ‘Doppler  bats’,  such  as  the  horseshoe  and 
pipistrelle,  while  chasing  insects  [9], 
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5.  Extensions  to  related  applications 

As  an  introduction  to  chase  problems,  this  paper  restricted  consideration  to  the  simple  case 
of  a  quarry  travelling  at  constant  velocity  followed  by  a  pursuer  travelling  at  constant  speed. 
In  this  concluding  section,  a  number  of  variations  on  this  basic  theme  are  listed.  Interested 
readers  are  referred  to  the  literature  for  further  discussion  of  them.  Many  of  them  could  form 
the  basis  for  challenging  student  projects,  well  suited  to  computer  simulation. 

Military  pursuit.  Anticipating  the  prey’s  future  location;  targeting  when  firing  bullets  at  the 
quarry;  limiting  the  maximum  acceleration  (minimum  radius  of  curvature);  pursuit  from  a 
constant  angle  of  attack;  optimum  search  trajectories  when  a  quarry’s  location  is  only  partially 
known;  controlling  the  impact  direction  so  that  it  is  from  the  front  or  side;  three-dimensional 
pursuit  through  the  air  or  water. 

Gaming  strategies.  Chase  in  a  confined  space  (such  as  a  boxing  ring);  pursuit  when  the  motion 
is  constrained  (such  as  on  a  chess  board);  preventing  a  quarry  from  reaching  a  refuge  (as  in 
tag);  a  runner  launching  at  a  moving  target  (as  in  dodgeball). 

Animal  motions.  Bugs  on  the  vertices  of  a  polygon  each  pursuing  their  nearest  neighbour  in 
a  clockwise  inward  spiral  [10];  a  dog  in  flowing  water  who  is  following  a  runner  on  land;  a 
spider  scurrying  along  his  web;  a  farmer  trying  to  round  up  multiple  escaped  pigs. 

Mathematical  analyses.  Polar  coordinates  of  the  prey  relative  to  the  pursuer;  quarry  following 
a  nonlinear  plane  curve  and/or  whose  speed  varies  in  a  predictable  manner;  relationship  to 
special  curves  such  as  the  tractrix  (shape  of  an  initially  straight  chain  whose  end  is  dragged 
in  a  direction  other  than  its  initial  axis)  or  caustic  (envelope  of  reflected  rays  from  a  curved 
mirror);  geodesic  pursuit  on  a  curved  surface  (e.g.  a  cylinder  or  hillside);  chase  on  a  rotating 
turntable. 
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